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Abstract 

The effective action for hard thermal loops in Quantum Chromodynamics (QCD) is 
related to a gauged Wess-Zumino-Novikov-Witten theory. Some of the technical issues of 
this approach are clarified and the Hamiltonian formulation is presented. The two-point 
correlator for the induced current in QCD is obtained; some simplifications of the dynamics 
of the longitudinal modes are also pointed out. 



I. Introduction 



The possibility of obtaining the quark-gluon plasma phase of hadronic matter in heavy 
ion collisions or perhaps in certain astrophysical contexts provides an adequate physical 
motivation for studying Quantum Chromodynamics (QCD) at high temperatures. In a 
perturbative analysis, hard thermal loops have to be calculated in order to carry out 
the resummation or rearrangement of thermal perturbation theory necessary to achieve 
a consistent expansion in powers of the coupling constant. Hard thermal loops are also 
part of the effective action for the long wavelength excitations of the plasma. They have 
therefore been the focus of many recent investigations [1-6]. 

The generating functional T[A] for the hard thermal loops is a gauge-invariant nonlocal 
functional of the gauge potential and is essentially an electric mass term for the gluons. 
The current given by T[A] leads to a gauge-invariant description of Debye screening and 
Landau damping effects [3]. r[^] also has an elegant mathematical structure, being closely 
related to the eikonal for a Chern-Simons theory [2] . The equations of motion, incorporat- 
ing contributions from T[A] are nonlocal, but it is possible to rewrite them as local, but 
coupled, equations [3]. Further, one can introduce an auxiliary field G, which has values 
in the color group, and an action r[74, G] which gives the local equations of motion [5]. 
The elimination of G leads back to the nonlocal functional T[A]. T[A, G] is related to the 
gauged Wess-Zumino-Novikov-Witten (WZNW) theory. A useful result which emerged in 
a rather simple way within this formulation is the positivity of the Hamiltonian in the hard 
thermal loop approximation. Actually, one can go further. The local nature of r[74, G] 
makes it a convenient starting point for the inclusion of hard thermal loop effects in the 
systematic calculation of any physical quantity. With this in mind, in this paper, we shall 
carry out the Hamiltonian analysis for this effective action. This will lead to a current 
algebraic description of hard thermal loops. A simple framework for the computation of 
current correlations will also emerge from this analysis. The two-point correlation function 
for the induced currents in the hard thermal loop approximation, which is related to the 
gluonic structure functions of the plasma, is obtained to the lowest order. We also point 
out some simplifications for the dynamics of the longitudinal excitations. 
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In section II, we give a brief resume of the effective action, set up the Hamiltonian 
analysis and obtain the operator algebra. In section III, the framework for computing the 
correlators and the simplifications for the longitudinal modes are discussed. There are 
three appendices; the first discusses the issue of whether the coefficient of the effective 
action should be quantized, the second with the definition of the Hamiltonian and the 
third with canonical derivation of the operator algebra. 

II. The Effective Action, Hamiltonian and Current Algebra 

The gauge potential = (— where t"' are hermitian N x AT-matrices repre- 
senting the generators of the Lie algebra of the color group; we take the color group to 
be SU{N). In addition to the gauge potential A^, we need the auxiliary field G{x,Q). 
This field is an S't/(A^) -matrix depending on spacetime coordinates and also on a unit 
three-dimensional vector Q. Further, G satisfies the condition G{x, —Q) — G^{x,Q). We 
shaU also use the notation Q^' = (1, Q), Q'^' = (1, -Q), and A+ = \A- Q, A- = \A- Q' . 
The effective action can be written as [5] 

S = j -\F^ + kV[A,G] (1) 

r[A,G] = j dVt j d^x'^S{G) + \ j d^xTi{G-^d-GA+ 

- A_d^GG-^ + A+G-^ A_G - A+AJ)] (2) 

where the integration indicated by dO, is over all orientations of the unit vector Q. S{G) 
is the WZNW-action for G, given by [7] 

S{G) = -— [ d^xTr{d+GG-^d-GG-^) + ^ [ TridGG'^ (3) 
27r Jm^ 127r J^a 

The integration in the first term is over the lightcone defined by the variables x'^ = 
x^Q^i x~ = x^Q'^; the second term, the so-called Wess-Zumino (WZ) term, requires 
extension of the fields to include one more coordinate, say s, and corresponding integration. 
We can take = x [0, 1] with fields at s = 1 corresponding to spacetime. Different 
ways of extending the fields to s 1 will give the same physical results. Since S{G) has 



integration over the lightcone variables, the integration for S{G) in equation (2) is over 
the transverse coordinates x'^ , ■ Q = 0. T[A,G] is gauge invariant and, if we ignore 
the dependence on the transverse coordinates, is just the standard gauged WZNW-action. 
Finally the constant k in equation (1) is {N + ^Np)^, where Np is the number of quark 
flavors and T is the temperature. The coupling constant g is not explicitly indicated; it 
can be recovered by the replacement — > gAfj^. 

The WZNW-action S{G) obeys the combination rule 

S{hG) = S{h) + S{G) - - [ Tr{h-^d_h d+G G-^) (4) 

Using this rule, the equations of motion can be easily derived as 

d+A- - d-A+ + [A+,A-] =0 (5) 
{D^F^'^'Y - J^" = (6) 

j-« = -A j dJ]TV{(-it«)[(^_-A_)Q^ + {A+-A+)Q"']} 

= _A / (mTi:{{-ie)[G-''D_G Q"" - D+G G-^Q""]] (7) 
ztt j 

D/j^ denotes the covariant derivative, D^'^ = d^5°''^ + f'^^'^A^^ and A+ = GA+G~^ — 
d+GG-\ A- = G-^A_G + G-^d_G. 

The current given by equation (7) agrees with what is obtained by evaluation of hard 
thermal loops, once the solution for A± from equation (5) is substituted in. The action 
(2) is thus an acceptable effective action at least at the level of the equations of motion. 
Actually one can make this a stronger statement. The equation of motion for the field 
G has no independent solutions. This can be seen as follows. We can parametrize the 
potentials A± in terms an SU {N)-m.atTix. V{x, Q) as 

A+ = -d+V V-\ A_ = -d-V V'-^ (8) 

where V'{x, Q) = V{x, —Q). The general solution to equation (5) is then given by 

G{x,Q) = V{x,-Q)B{x+,x'^,Q)C{x-,x'^,Q)V-\x,Q) (9) 
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where C is an arbitrary 5'L/"(A?")-matrix depending on the variables indicated and B is 
given by C with Q ^ —Q. The matrices B, C represent the new or independent degrees of 
freedom for the field G. Notice, however, that the paramctrization (8) of the potentials has 
redundant variables; a transformation V VU{x~ , x'^ ,Q), with a corresponding change 
in V', leaves the potentials invariant. Further, ordinary gauge transformations act on the 
matrix F as F — > h{x)V. Thus the physical subspace for the components A± is given in 
terms of the matrices V with the identifications 

V{x, Q) ~ h{x)V{x, Q)U{x-,x'^, Q) (10) 

The gauge freedom of multiplying V^s by matrices which do not depend on x~^, viz. C/'s 
in (10), shows that the we can reduce G to just V{x, Q)V~^{x, Q). There are thus no new 
real dynamical degrees of freedom in G. The action can be simplified for this G to 

T[A,G] ^ - J dn d^x^ S{V-\x,-Q)V{x,Q)) (11) 

This agrees with the expression obtained by Feynman graph evaluation of the hard thermal 
loops. 

The Hamiltonian for the effective action (2) was given in ref.[5] as 

n= jd^x I ^ ^ ^/ d^^li^oG DoG-^) + {Q-DGQ-DG-^] -^[Jg^^j 

(12) 

where = E^, F«. = e^-fc^" and 

= (D-E)" + ^ f dnTT[{-it''){G-^D_G - D+GG-^)] (13) 

= is the Gauss law of the theory; it is also the time component of the equation 
of motion for the gauge field. Expression (12) makes it clear that the Hamiltonian is 
positive for all configurations which are physical, i.e. obey the Gauss law. (It may be 
worth recalling that, even for the Maxwell theory, the canonical Hamiltonian is positive 
only for physical configurations.) 

The Hamiltonian analysis is most easily carried out, as in the usual cases, in the gauge 
= 0. We start by defining the currents 

J+ = ^D+GG-^ = (-it'')Jl 
47r V / 



J_ = -^G-^D_G = {-1^)31 (14) 
47r 

By virtue of the property G~^{x,Q) = G{x,—Q), these are related by J+{x,—Q) = 
J-{x,Q). The Hamiltonian can be written in terms of these currents as 



(15) 



We have chosen the Aq = gauge; Gauss law must henceforth be imposed as a constraint. 
( For fixed Q, the integrand of the second term involving the square of J+ and the square 
of J_ is the Sugawara form of the Hamiltonian, well-known for two-dimensional current 
algebras.) The equations of motion in the Aq — gauge are 

Ef =doA^ (16a) 
doEf + eijk{DjBkr = j dil Qi{Jl - Jl) (166) 

{D-J+T = -^^tQi (16c) 

Equation (16a) is just the definition of the electric field; however, in a Hamiltonian ap- 
proach it is an equation of motion and we have displayed it as such. 

The commutation rules must be such that equations (16) follow as the Heisenberg 
equations of motion for the Hamiltonian (15). Knowing the current algebra of the WZNW- 
theory, we can make a guess as to what the appropriate commutation rules are for our 
problem and verify them by checking that they lead to equations (16) starting from the 
Hamiltonian (15). The commutation rules are then seen to be 

[ E'l{x) , A]{x') ] = -iS^^SijSix - x') (17a) 

[Et{x),Ji{x')] =±i^Qi5''''5{x-x') (176) 

47r 

[ J5(x,Q),4(f' ,g') ] =ir^'^Jl5{x-x')5{Q,Q') 

T^Qi{D,)fS{x-x')5{Q,Q') (17c) 
[Jl{x,Q),j'_{x',Q')] =0 {17d) 

All other commutators vanish. S{Q, Q') stands for the 5- function on the two-dimensional 
sphere corresponding to the unit vector Q, i.e. J dCl^, 5{Q^Q')f{Q') = f{Q). 
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We have checked that these rules obey the Jacobi identity. Since we are postulating 
the algebra, this is a necessary check. Of course, commutation rules can be obtained 
from the action by standard quantization procedures. We have verified that equations 
(17) follow from such a canonical approach as well. The canonical method is somewhat 
involved. The most eflicient way is to use the language of symplectic structures. This is 
sketched in the appendix. Finally, the condition J^{x — Q) — J^{x,Q) = has to be 
imposed as a constraint, just like the Gauss law. 

We now turn to deriving an equation for the correlators of currents. The parameter 
k in equation (17c) controls the scmiclassical expansion, large k being the classical limit 
[7]. For a scmiclassical expansion, it is convenient to use the currents I± = ^J~^J±. 
Eventhough we are not doing a scmiclassical expansion at this stage, it is useful to consider 
correlators of J+'s. (7_'s are related by Q — > Q')- Thus introducing a source a°'{x, Q) for 
/^(x, Q) and a source rj^{x) for A^^x), we can define a generating functional for correlators 

Z[r7,a]=exp(^^[r/,a]) = (Texp[^ j d^x {{ J dil a^I^) + r^fA^}]) (18) 

^[77, q;] generates the connected functions. The angular brackets denote a thermal average 
or a vacuum expectation value depending on whether the thermal effects of the soft modes 
(momenta ~ gT) are to be included or not. We now derive Z[r], a] with respect to a and 
take the time-derivative. As is well known, this leads to a term involving the commutators 
and a term involving the time-derivative of the current I^. This can be simplified using 
the equations of motion and the operator algebra (17). The result is 

{^o + Q^V,){^{x,Q))+ r^'^(^^a\r) + g,(A,^(x)F(a;,Q)) 

+ J ^Qido{A1) - {QiViaY =0 (19) 
V 47r 



where 



(/"(x,g)) = (-i / )w, 

S . , . S 



(-4»Wr(..Q)) = H^)(-i^^)H.. etc. (20) 
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The sources are not yet set to zero in equations (19,20). Further, Df'^ = di5°''^ + /"^'^(A^). 
In a similar way, we can derive an equation for the quantity {A'^[x)). This is given by 



where 



Vt =0 (21) 



/ 



abc 



dj{A'l{x)A'^j{x)) + {A){d,A^^ - d^A^)) + r\{A){x)Al{x)){A^,{x)) 

+ {A]{x)A]{x)){A\{x)) + {A]{x)){A\{x)A]{x)) 



+ 



{A){x)A\{x)A^^{x))) 



(22) 



and J^. = di{A^) - d^{A1) + f^'^ {A\) {A<^) . 

Any correlation function involving the fields can be calculated using the equations 
(19-22). It is, of course, possible to write functional integrals for the correlation func- 
tions. However, since the currents of the WZNW-theory are the important quantities for 
calculating the correlators of the fields A^^ the present approach of obtaining differential 
equations using the Hamiltonian framework is somewhat simpler. Equations (19-22) are 
written without the renormalization factors; eventually, each average of an independent 
opertor product should have such a factor. 

As an application of these equations, we calculate the correlator {I"'{x,Q)I^{y,Q')) 
to the lowest order (with no loops). Differentiating the equations (19,21) with respect to 
<^^{yjQ') and setting all sources to zero, (indicated by a subscript zero on the angular 
average signs), we get 



[(^2 - dkdk)d,j + d,dj] {A^ix)l\y, Q'))o 

iQ;rt(x-y)-2 J dn" Q':{I%x,Q")l\y,Q'))o 



An 



C(23) 



{do +Qid,){r{x,Q)l\y,Q'))o + ^ ^do{At{x)l\y,Q'))o 

+ i{Qidi)5''^5{x - y)5{Q, Q') ^ 0. 



(24) 



We have dropped the higher loop terms, as emphasized by the approximate equahty. Also 
we note that it is consistent to take = 0. Equation (23) can be solved by using the 
Green's function 



-I 



Po 



(25) 



The integrand is singular as po ^ 0, indicating the residual gauge invariance under time- 
independent transformations. Nevertheless we can use this in the usual way with some 
prescription for handling the singularity such as the a-prescription [9] . Using the solution 
of (23) in (24) we get the equation 



/ 



{do + Qidi)S{x - z)S{Q, Q") + ^do{QiGij{x, z)Q'^) 



= S 



{I%z,Q")l\y,Q')) 
(26) 



-iQ^d,d{x - y)S{Q, Q') + —^o{Q^Gij{x, y)Q'^) 
The solution to this equation is given by 

(r(x, Q)l\y, Q')) = I ^^e-'^^-y) F{p, Q, Q') (27) 
F{p. g, Q') = h5{Q, Q') +h\Q- Q' - ^' ^' ] + f /' ^ t ^' (28) 



Pi 



Po 



fi 



ip- Q 



{po-p- Q) 
1 



_ Po+p-Q 



f - 2^ ^ 
Po 



{po-p-Q-la (g)) 



Att p"^ 



(29) 



We have used the vacuum average. A thermal distribution for the soft modes can be taken 
into account by imposing periodicity in imaginary time. The correlator (27) is useful, for 
instance, in calculating the scattering of a gluon in the quark-gluon plasma. It is what 
defines the gluonic structure functions of the plasma. 
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We now turn to a slightly different topic. The Hamiltonian analysis and the equations 
of motion can be considerably simplified for the longitudinal modes. Eventhough the re- 
striction to just the longitudinal modes is not physically well motivated, the simplifications 
are interesting and could eventually be incorporated in a more general framework. For this 
reason, we shall briefiy discuss the reduction to longitudinal modes. 

The longitudinal modes can be parametrized in terms of a matrix field U as 

Ai = -diU U-^ (30) 

U is in general time-dependent and hence the gauge potential (30) is not a pure gauge. 
The magnetic field, however, is zero. The electric field can be parametrized in terms of a 
field n as £'i = —U{diIl)U~^ . We can also define currents 

3± = U-^J±U (31) 

j± and n are gauge invariant degrees of freedom; the field 11 is constrained in terms of j± 
by the Gauss law. Since J± commute with Ai, they obviously commute with U. Using 
this fact, the operator algebra can be simplified. The nontrivial commutators become 

[m^''ix),djdj-n\y)] =ir^'dkdkn%x)S{x-^ (32a) 

[ -d,d,U-{x),jiiy,Q) ] = -irf"^jl[x,Q)Six-y) 

±^5'"'QidiS{x-y) (326) 

[m^,Q),ji{y,Q') ] =ir'''mx,Q)6{x-mQ,Q') 

T^S'^'QidiSix - y)S{Q. Q') (32c) 

The Hamiltonian becomes 



j (fx \{diWdiW) + dn{jXj% + 



(33) 



The Hamiltonian does not involve the field U directly and hence the commutator of H" 
with U, which is nontrivial, is not needed to analyze the dynamics of the longitudinal 
modes. 
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The Gauss law condition becomes 

-dATl'' + JdQ iJl + r_) = (34) 

The commutation rules (32) are consistent with imposing the Gauss law as a strong 
operator condition. We may thus take (34) as the definition of and the Hamilto- 
nian and the commutation rules (32c) define the dynamics, with one constraint, viz. 
jl{x-Q)-r-{x,Q) = 0. 

The Hamiltonian equations of motion are now 

doC^ - QidiCS + r^'C^oCI = (35) 

where = ^^i^ Cq =Tl°- — C^. This equation can also be written as 

{do + Qidi)Ct + r^'^n^Cf = QidiW' (36) 

We have a similar equation for C'^ = ^jl, obtained from (36) by Q — > —Q- The Gauss 
law is 

-dAU^ + ^Jdn {CI + ) = (37) 

The solution of these equations, if we neglect the term involving the structure constants 
f°'^'^, are the Abelian longitudinal plasma waves. Equation (36) can thus be used to 
systematically improve on the Abelian approximation. 



Appendix A 

As is well known, WZNW-theories for unitary matrices lead to a quantization require- 
ment on the constant k multiplying the action. This arises from the single-valuedness of 
the transition amplitudes or wavefunctions. For our case, eventhough the auxiliary field 
G is a unitary matrix, there is no such requirement. The Wess-Zumino term which leads 
to the quantization requirement can be written, for a fixed choice of Qi, as 

Swz = ^ f d^x^ I Tr(dG Q-^f (Al) 
127r J Jm3 
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where is a space whose boundary is the two-dimensional world of the lightcone co- 
ordinates x^. Swz can be thought of as being obtained by integrating 5Swz with the 
extra coordinate of M^, say s, parametrizing the path of integration. Physical results are 
independent of this path. The difference for two extensions of the field G into the extra 
dimension corresponding to s is given by 

This is for the intrinsic two-dimensional problem, ignoring the transverse coordinates. 
is the three-dimensional sphere. Q is an integer corresponding to the winding number of 
the map G over the three-sphere. (The nontriviality of (A2) has to do with the fact that 
Tr(dG G~^)^ is a closed but nonexact differential form.) The single- valuedness of e*'^'^^^ 
then requires k to be an integer. 

For our case, we have the two transverse coordinates as well; we can write equation 
(Al) as 

Swz = tI- [ T^idG G-y CO (A3) 

where u = eij^Q^ dx^ dx'^ . There is no restriction for the exterior derivatives in (A3) to 
be with respect to x'^ or s; this is enforced by the choice of to. Notice that eventhough 
Tr(ciG G~^)^ is not exact, the differential form in (A3) will be, if to is exact, uj can be 
formally written as dp, but the one-form p is linear in the coordinates and so we cannot use 
this unless we impose some stringent fall-off conditions on the fields G. However, when we 
take the difference of two extensions of G into M^, we get the quantity on the right hand 
side of (A3) integrated over a five-sphere, p exists with no singularity and the integrand 
is exact. There is no nonzero winding number and no quantization of k. 

Appendix B 

The Hamiltonian for the effective action was obtained in ref.[5] by taking the derivative 
of the action with respect to the final value of time in a discretized version. This is a 
standard result, but does not seem to be well appreciated. A short argument justifying 
the calculation will be sketched here. 
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We shall use the example of a scalar field to illustrate the calculation. Consider the 
Lagrangian 



dt 



where includes terms which do not depend on time-derivatives of 0. The Hamiltonian, 
by definition, is the generator of reparametrizations of the time-variable t. Let us consider 
therefore the infinitesimal transformation t ^ t = t -\- e(t), where e{t) is a function of t. 
We then have di = dt{l + e'), = {1 — e')dt, where the prime denotes differentiation with 
respect to t. For the variation of the action S = J dt L, we then find 



5S 



d'^x 



d(f) 
'di 



+ 



E S(f) 



5(f) is the change induced in and £ is the variational derivative. (Setting £ to zero will 
give the equations of motion.) The final limit of the time- integration is not changed in this 
variation, only the coordinate label for it is varied. From equation (B2) we can identify 
the Hamiltonian as 



d^x 



dt 



+ v{4>) 



m 



If one imposes the equation of motion, viz. |^ = 0, one can prove the conservation of H. 
However, to identify the expression for H we do not impose the equations of motion; we 
just use the separation of 5S in the form given in (B2). 

We can recast the above result as follows. Let us write the integral over time by 
discretizing the time-interval, with tN being the final time-label. Then e'itn) = , 
etc. The action can be written as 



S= d-^x 



(A 



■'n—i — lj 



{tfi—i tn—i—l) 

The expression (B2) for 5S can be written as 



{V {<!>)) ^_,{tn-i-tn-i-l) 



m 



5S = - [(ejv - eN-i)H{tM) + ...] + j£ 



{B5) 



We see that H can also be identified as the coefficient of ejv in the first set of terms. 
However, sn is the shift in tjv; thus the coefficient of ejv is the derivative of the action with 
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respect to the final time-label t^. However, in taking this derivative we should not vary 
the fields as functions of time so as to avoid contributions from the second set of terms, viz. 
those proportional to £. This qualification means that we can use expression (B4) for the 
action, but interpret (pn-i as independent variables, not as (pitn-i). Further, although we 
change tjv, the number of segments of time interval is fixed. The limits of the summation 
are unchanged, reflecting the fact that the variation in (B2) only changes the coordinate 
labels. Applied to the effective action, the above argument gives the Hamiltonian (12), i.e. 
the result in ref.[5]. 

Appendix C 

In this appendix, we discuss the canonical derivation of the operator algebra. 
The general method is as follows. Given an action S depending on a set of fields 
</7r, r = 1, 2, ...n, we consider its variation under a change of fields. This is of the form 



where £r is the variational derivative. The second term is the surface contribution, with 
K^^ being linear in the variations 5(pr- If we consider the time-coordinate as the evolution 

parameter and impose suitable boundary conditions at spatial boundary (viz. that Ki 
vanishes there), we can write the surface term as the difference of ©'s at the initial and 
final time-slices, where 



is the canonical one-form. The variation of 0, antisymmetrized in the variations, (i.e. 
the exterior derivative on the field space ) gives the symplectic structure to. (This is 
an old result going back to nineteenth century analytical mechanics; in the context of 
quantum field theory, Schwinger's action principle is closely related to this. For some 
modern references, see [10]. Also, an analogous method can be used even if the initial and 
final data hypersurfaces do not correspond to constant tinic-sliccs.) There arc ambiguities 
in identifying © as above, since we can add total divergences to the Lagrangian; ui has 
no such ambiguity. (The changes in © given by such total divergences are, of course, the 
canonical transformations.) 




{CI) 




(C2) 
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Given cu, the Poisson brackets for the canonical coordinates are given by the inverse 
to coij, viewed as an antisymmetric matrix, where lu = ^Uijda'^ 5cr^ , and cr* denote the 
canonical variables (like the fields or their time-derivatives). The commutators are then i 
times the Poisson brackets. For the present problem, we need commutators for currents, 
rather than the canonical coordinates, and for this purpose it is algebraically simpler 
to use an equivalent definition. Given u, a vector field is said to be Hamiltonian if 
Ca\^ = —SA, for some function A; alternatively, given any function A, we can associate a 
vector field to it by this condition. The notation is as follows. If we write u in terms of 
its antisymmetric components ujij and ^ — C jfi > then 

^Ja; = fa;y5(7^' (C3) 

The commutator of two functions A and B is given by 

[A,B] =iiU\^B\oj) (C4) 

This method can be applied to our problem. The only technical complication is that 
the WZ-term requires an extension of the fields into an extra dimension. We can take the 

extended space as x [0, 1], s = 1 being real spacetime. The variations of the fields 
at s = can be taken to be zero. From the variation of the effective action and after 
identifying 0, the symplectic structure can be obtained as 

u; = J (fx 5Ef5A1 + ^ J dn Tr[v SLq + v^{Lo + QiLi) 

+ vQi5Ai + G-'^vG QiSAi] (C5) 

where v = 5G G~^ , = d^G G~^ , = G~^d^G. The Hamiltonian vector fields are 



For the variation 6G — 9G, the vector field is 



= V{e) + J (fx dQ ^Qie\5''' - 2M^«)^ {C7) 
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where M«'' = TT[{-it'')G{-it^)G-^] and V{d) replaces v hj 9 and 6Lo by [9,Lo + QiU]. 
In other words, 

We can calculate the various commutators using these vector fields and the formula (C4). 
The result is the algebra (17) in the text. 
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